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The mystery of the normal state in the underdoped cuprates has deepened with the use of newer 
and complementary experimental probes. While photoemission studies have revealed solely 'Fermi 
arcs' centered on nodal points in the Brillouin zone at which holes aggregate upon doping, more re- 
cent quantum oscillation experiments have been interpreted in terms of an ambipolar Fermi surface, 
that includes sections containing electron carriers located at the antinodal region. To address the 
question of whether an ambipolar Fermi surface truly exists, here we utilize measurements of the 
second harmonic quantum oscillations, which reveal that the amplitude of these oscillations arises 
mainly from oscillations in the chemical potential, providing crucial information on the nature of 
the Fermi surface in underdoped YBa2Cu30e+x. In particular, the detailed relationship between 
the second harmonic amplitude and the fundamental amplitude of the quantum oscillations leads 
us to the conclusion that there exists only a single underlying quasi-two dimensional Fermi surface 
pocket giving rise to the multiple frequency components observed via the effects of warping, bilayer 
splitting and magnetic breakdown. A range of studies suggest that the pocket is most likely associ- 
ated with states near the nodal region of the Brillouin zone of underdoped YBa2Cu306+ x at high 
magnetic fields. 



While high-T c cuprate superconductors have a history 
spanning a quarter of a century 0, Q , quantum oscilla- 
tions in these materials have only been measured as re- 
cently as a few years ago Interestingly, their advent 
has signalled a re-evaluation of the electronic structure in 
the normal state induced by high magnetic fields in the 
underdoped cuprates [il-fl2j]. Where previously, the very 
existence of quasiparticles was questioned in the normal 
state of the underdoped cuprates [liij . quantum oscilla- 
tions have uncovered quasiparticles consistent with the 
existence of a Fermi liquid within this regime Q . 

Despite the progress made by quantum oscillations 
that complement a broad array of other experiments in 
discerning the normal state of the cuprates, fundamental 
mysteries remain both in interpreting the results of quan- 
tum oscillation measurements, and in reconciling them 
with other experimental techniques. Quantum oscilla- 
tion results have been interpreted in terms of an ambipo- 
lar Fermi surface with multiple pockets @, H, H, E, 12 1 
(i.e.) comprising hole carriers concentrated at the nodes 
] and electron carriers located near the 
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antinodes k an ti nod ai = [f ,0] & [0, f], where a is the lat- 
tice constant. This picture, however, is in direct contrast 
with the observation made by photoemission experiments 
of a significant (~ 50 meV) gap in density of states at 



the antinodes 14-1 



An unambiguous interpretation of quantum oscillation 
results in terms of charge carrier type, and single or mul- 
tiple surfaces in the case of the underdoped cuprates has 
thus far been prone to pitfalls, due to the quasi- two- 
dimensionality of the Fermi surface, closeness of frequen- 



cies and effective cyclotron masses of the low frequency 
components, and the lack of clarity regarding the ob- 
servation of a significantly larger frequency & E3- 
Quantum oscillations in the case of a fixed chemical 
potential cannot discern the sign of charge carriers or 
whether there is a single or multiple types of charge car- 
rier present [lj|, causing recourse thus far to less di- 
rect probes such as the Hall effect and thermopower, 
and calculations of the electronic structure to infer a 
Fermi surface comprising mobile electrons and less mo- 
bile holes jiSiJii. 

In this work, we observe well-defined second harmonic 
quantum oscillations from a broad array of torque and 
contactless resistivity experiments performed in a 45T 
DC magnet, a 60T generator-driven magnet, and 65 & 
85T pulsed magnets in angles J$ # J$ 60° and tempera- 
tures IK < T < 5 K in the National High Magnetic Field 
Laboratory. From the analysis of the measured second 
harmonic oscillations, we present the surprising finding 
that the harmonic content in underdoped YL^CusOg+x 
(x = 0.56) arises from oscillations of the chemical po- 
tential, enhanced compared to that in normal metals 
by t he q uasi-two-dimensional topology of the Fermi sur- 
face [3 l23l426| . Oscillatory features corresponding to 
multiple frequencies in 1/B previously reported 0,SI3i0| 
are observed: (F a = 535(5)T, F yl = 440(10)T, F l2 = 
620(20)T) down to 22T, and F f3 = 1650(20)T at the low- 
est temperatures [27| : yet we conclude from the magni- 
tude and phase of the second harmonic of the chemical 
potential with respect to the fundamental oscillations, 
that the Fermi surface consists only of a single quasi-two- 
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dimensional pocket, and no other Fermi surface section of 
significant thermodynamic mass, the multiple frequency 
components arising from effects of finite c-axis disper- 
sion, bilayer splitting, and magnetic breakdown. Given 
the nodal concentration of the density of states observed 
in zero field photoemission (l^ - [l6| and high field heat 
capacity experiments 17|, indications are that, in the 
absence of a significant effect of magnetic field, the ob- 
served Fermi surface pocket is located at the nodal region 
of the Brillouin zone. This finding does not, however, 
necessarily imply a contradiction with the observation of 
a negative Hall and Seebeck coefficient 0, [27| ■ 



RESULTS 

A number of observations relating to the second har- 
monic reveal departures from the oscillatory behavior ex- 
pected from the standard Lifshitz-Kosevich (LK) theory 
that assumes a fixed chemical potential [3] . Figure [I]a 
shows the second harmonic extracted from quantum os- 
cillations measured on YBa 2 Cu306+ x at a series of fixed 
angles 9 (between the crystalline c-axis and the magnetic 
field) as a function of the c-axis magnetic field compo- 
nent Bcos9, where B = pioH is the applied magnetic 
field (up to 65T in a pulsed field magnet). We find that 
the second harmonic oscillations retain a constant phase 
between 9 = 0° and 57° at fixed B cos 9 A remarkable 
feature in the data is the absence of a zero near w 18°, 
which is expected in the second harmonic if fj, is fixed, 
given the presence of a spin zero in the fundamental at 
« 54° [27J (see Fig. [1}d). Instead, the second harmonic 
amplitude falls to zero at precisely the same angle as the 
fundamental. 

Figure [Tt compares the T dependence of the funda- 
mental amplitude A\, with that of the second harmonic 
amplitude Ai on a linear-log plot normally used to infer 
the appropriate cyclotron mass. These precisely temper- 
ature controlled measurements were performed with the 
sample in liquid 4 He in a 60T generator-driven magnet. 
Quantum oscillations corresponding to a fixed chemical 
potential are expected to exhibit a second harmonic with 
a cyclotron mass twice the fundamental cyclotron mass, 
yet strangely a fit of the second harmonic to the thermal 
damping factor (see Fig. QJ:) yields a mass that is too low 
by about 25 %. 

Figure [2] shows the detected amplitude of the second 
harmonic as a function of magnetic field measured using 
two different methods: magnetic torque in DC fields and 
contactless resistivity in pulsed fields. The measured sec- 
ond harmonic amplitude is found to scale precisely as the 
square of the fundamental as a function of magnetic field. 
This square relationship mirrors that as a function of an- 
gle (Fig.Q]D) and as a function of temperature (Fig.QJ:). 



observed second harmonic oscillations and the fundamen- 
tal oscillations at 3.8° in the measured torque, 0° in the 
contactless resistivity signal and at 38° in the contactless 
resistivity. The second harmonic retains a constant phase 
with respect to the fundamental, consistent with the saw- 
tooth waveform shown schematically in Figs. [3J1 and e 
respectively. This is inconsistent with that expected for 
quantum oscillations corresponding to a fixed chemical 
potential, which leads to a different form as shown in 
Fig. EJf) and to a spin zero and phase inversion at 18°, 
which we do not observe. 

Finally, Fig. @Ji shows a collation of quantum oscilla- 
tions measured in magnetic torque and contactless resis- 
tivity, spanning overall a range 22T to 65T. The identi- 
cal beat structure previously revealed to correspond to 
three closely spaced frequencies (F a = 535 (5)T, F 7 i = 
440(10)T, F 72 = 620(20)T) is found to persist down to 
22T, confirming the multiple frequency components ear- 
lier reported [1 @, 0, @, El • 



DISCUSSION 

We broadly consider three chief aspects of the exper- 
imental data, and understand the unusual behaviour of 
the second harmonic in terms of the effect of oscillations 
of the chemical potential [i. 

Large oscillations in /i occur in layered metals owing to 
the absence of a large reservoir of phase-incoherent states 
to pin /i at an essentially constant value [231 - 125] . Instead 
of the Landau levels passing through a constant chemical 
potential — as described in the simplest LK theory [19( — 
fi itself oscillates. The chemical potential remains pinned 
to the highest partially occupied Landau level as its de- 
generacy changes, jumping from one to the next as the 
levels become sequentially depopulated on increasing the 
magnetic field (shown schematically in Fig.[3}i). Oscilla- 
tions in the magnetization (the de Haas-van Alphen ef- 
fect) are affected by the jumps in /i, yielding a waveform 
that resembles a sawtooth (see schematic in Fig. [3^). 

We consider the effect that the oscillatory component 
jl has on modulating the Onsager phase of the quantum 
oscillations <p = 2nF/B — rj (where r\ ss tt) arising from 
a single two-dimensional pocket giving rise to quantum 
oscillations of frequency F . The oscillatory magne- 
tization is given by 
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Fig. [3] allows us to infer the phase relation between the 



where n p is the product of quantum oscillation ampli- 
tude damping factors (e.g. Dingle, thermal, spin, mosaic 
spread, etc.) affecting the p th harmonic [19(, and Mq 
is the absolute amplitude at infinite field. The effect of 
chemical potential oscillations are explicitly included in 
leading order by the expansion of <fi([i) = tfio+fiidtfto/dfl), 
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FIG. 1: a, Field-dependence of the second harmonic (2F) from the data in Ref. [271 ] at different angles 9 between the crystalline 
c-axis and the magnetic field after having subtracted fits to the fundamental (F) and a slowly varying background. The data 
are plotted versus BcosO to facilitate a comparison of the relative phase. The second harmonic amplitude falls to a small 
value near 54° (at which the fundamental undergoes a spin zero [27|), following which it recovers with the same phase at 
higher angles, b, A comparison of the Fourier amplitude of the second harmonic obtained over the field interval indicated 
against the rescaled square of the amplitude of the fundamental (F 2 ) reveals that they precisely track each other over a broad 
range of angles, spanning the spin zero in the fundamental near 54°. c, Temperature-dependence of the fundamental (black 
squares) and second harmonic (colored squares). We plot (1 /n) ln[A n B / £T] versus £T/B, where n is the harmonic index, 
£ = 2-K 2 k^m e /he = 14.69T/K and A n is the measured harmonic amplitude. The slope of this relation is expected to have the 
tti 

is the free electron mass, for both the second harmonic and the fundamental (given the division by n 
. We find, however, an unexpectedly low slope of value 1.3(1) for the second harmonic, which is lower than 
1.6(1) measured for the fundamental. On plotting (1/n) \n[A n (B / £T) n ] versus £,T/B for the second harmonic 



value , where m. 

on the left axis 
m* 



the slope 

(pink symbols, right axis), however, the slope of the relation agrees with that of the fundamental. The latter implies that the 
temperature dependence of the second harmonic is given by (where R T i = X/sinh X w 2Xe~ x for X = £,Tm*/Bm c < 1), 
as expected from the effects of oscillations in the chemical potential [28j |. 



where tf>o — 4>(p), and /i is the steady part of /x 28. 2 
This allows us to rewrite Eqn. [T] in the form 
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where T is determined by the relative strength of p, as 
discussed below (0, Hij]). While the fundamental fre- 
quency continues to be well described by the standard 
theory for oscillations corresponding to fixed /i, the sec- 
ond harmonic can be dominated by a new contribution 



of relative weight 



2Tni 

K2 



(see Eqn. [2]) resulting from the 



oscillations in /i. Here V represents the fraction of the to- 
tal density of states associated with the postulated two- 
dimensional pocket. 

In Fig. [4ji we compare the fundamental and sec- 
ond harmonic magnetization inferred from the measured 
torque M oc f/B and from the contactless resistivity 
M oc J pdB/B 2 [lj|, where the integration is over the 
experimental field range. Their similar behavior and seal- 
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FIG. 2: a, A comparison of the second harmonic oscillations (thick purple lines) with the scaled square of the amplitude of the 
fundamental (dots) versus magnetic field, indicating that they are closely correlated for the oscillatory magnetization M oc f/B 
extracted from the magnetic torque f at an angle 8 « 3.8° measured in the 45T DC magnet. The dots are extracted from the 
scaled square of the fundamental oscillations maxima and minima. The actual fundamental oscillations are depicted in light 
grey, b, The same repeated for the oscillatory contactless resistivity data at 8 ~ measured in the 65T pulsed magnet, showing 
the same behavior as in a. c, A comparison of the second harmonic oscillations (thick navy lines) with the scaled difference 
between the squared amplitude of the fundamental a and 7 components (a and 7 component fit from ref. 0|). This would be 
the expectation for a and 7 corresponding to opposite carrier types, we find that the observed second harmonic does not follow 
this form. 



ability indicates that the oscillations in both quantities 
have a closely related origin [bij . 

First we discuss the phase of the observed second har- 
monic oscillations. We see that in the case of strong 
chemical potential oscillations (where the first term in 
the prefactor of sin20o in Eqn. [2] is dominant), the sec- 
ond harmonic oscillations are in fact out of phase by 7r 
from the conventional case of oscillations corresponding 
to a fixed chemical potential — precisely as we observe in 
experiment (predictions and comparison with experiment 
shown in Fig. |3ji.) 

Secondly, we discuss the amplitude A2 of the second 
harmonic, which is expected to vary as the square of the 
amplitude A\ of the fundamental in the limit 2Yk\ > n 2 
in the case where harmonic enrichment arises mainly 
from chemical potential oscillations. The experimentally 
observed square relation A 2 oc A\ as a function of mag- 
netic field at various angles both in the torque measure- 
ments and contactless conductivity (Fig. [2] and Fig. Hp) 
is consistent with this expectation. Also consistent with 
this expectation is the observed absence of a low-angle 
spin zero as discussed earlier for the second harmonic 
(Fig. [T|d( and the low value of the effective mass inferred 
from a conventional linear-log plot of the second har- 
monic as also mentioned earlier (Fig. [I};) . 

Thus far we have considered a model with a single two- 
dimensional Fermi surface pocket and hence quantum 
oscillations with a single periodicity in 1/B. The pre- 
dictions of a model including multiple frequencies arising 
from warping and bilayer splitting of the same pocket are 
not substantially different from those discussed above. 



However, the inclusion of a second two-dimensional Fermi 
surface pocket associated with carriers of opposite effec- 
tive charge would lead to an additional harmonic enrich- 
ment term from chemical potential oscillations of the op- 
posite sign (i.e.) — nf in addition to — n\ in Eqn. [21 lead- 
ing to behaviour in general inconsistent with harmonic 
components scaling as the square of fundamental com- 
ponents [19} . Fig. [5J; shows the predicted behaviour for 
the second harmonic were the a and 7 component of the 
oscillations to be associated with pockets of opposite car- 
rier type — this does not fit the experimentally observed 
form of the second harmonic, ruling out the possibility of 
the observed quantum oscillations arising from pockets 
containing opposite carriers. 

Our finding A 2 oc A\ therefore supports the idea that 
the observed frequency components are all closely related 
and are likely to arise from the same two-dimensional 
pocket, its purely two-dimensional structure modulated 
by weak warping and bilayer splitting. We next con- 
sider what can be learned from the magnitude of the 
ratio of A2 to A\ , namely, from the magnitude of T as 
defined in Eqn. [5J This ratio is most accurately deter- 
mined at 9 = 38° where the 7 oscillatory component is 
weak and the quantum oscillations are dominated by a 
single periodicity due to the a component 

DOS(a) 



In this case T 



where 



DOS(a) + DOS( 7 ) + DOS(res) ' 
'DOS' stands for density of states and 'res' for any reser- 
voir of carriers at the chemical potential not associated 
with the a or 7 components. From the measured value 
r(38°) = 0.50(15) at this angle (see Fig. gp), we esti- 
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FIG. 3: a - c, The fundamental (green) and second harmonic (magenta) extracted from the magnetic torque measured in 
YBa2Cu306+a: at different field orientations 9 as indicated, revealing the same phase behavior as seen in the schematic figure 
(e). The magnetization is inferred using M oc f/B and M oc J pdB/B 2 from the magnetic torque f and contactless resistivity 
respectively. Integer Landau level filling factors (when fi resides in the gap between two-dimensional a component Landau 
levels) are indicated, d, A schematic showing the Landau level energies (cyan lines) and how the chemical potential /x (red 
line) and consequently the magnetization (red line in e) oscillates for a single two-dimensional electron pocket, /i jumps 
between Landau levels (cyan lines) at integer filling factors causing it to be momentarily situated in the Landau gap. e. The 
magnetization (red) together with the extracted fundamental (black) and second harmonic (magenta), f, The corresponding 
magnetization waveform (blue) when the chemical potential is fixed, giving rise to a second harmonic of opposite sign. Vertical 
lines in gray correspond to integral Laundau levels, and in orange correspond to half-integral Landau levels. 



mate an upper bound for 



DOS(res) 



of about 



DOS(a) + DOS( 7 ) 
30%. This suggests that our model of a single quasi- two- 
dimensional Fermi surface pocket, giving rise to multiple 
frequency components a and 7 together with possible 
magnetic breakdown components, provide an essentially 
complete description of the Fermi surface under our ex- 
perimental conditions in underdoped YBa2Cu306+x- 



Complementary photoemission [l44l6j and high field 
heat capacity [17| experiments yield a concentration of 
density of states at the nodal location in momentum 
space, suggesting that in the absence of a significant ef- 
fect of magnetic field, the single Fermi surface pocket that 
gives rise to multiple frequency components observed 



by quantum oscillations in underdoped YBa2Cu 3 0e+ :I ; 
is also connected with the nodal region of the Brillouin 
zone. 



CONCLUSION 

In summary, we find that second harmonic quantum 
oscillations in underdoped YBa2Cu306+ x arise mainly 
from the effect of oscillations of the chemical potential. 
An analysis of the phase and magnitude of the second 
harmonic components relative to the fundamental com- 
ponents allows us to conclude that the Fermi surface con- 
sists of a single quasi-two-dimensional pocket, the mul- 
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FIG. 4: a, Comparison of the magnetization M oc f/B and M oc f pdB/B 2 inferred from the measured magnetic torque f 
(purple line) and contactless resistivity p (red dotted line) respectively, b, Plot of the second harmonic (at 9 « 36.3° in fields to 
65 T and at 6 ~ 38° in fields 85 T) divided by the infinite field intercept po determined in (c): (fo/po) represented by solid lines 
is compared with (pi/po) 2 , represented by solid dots. Here the amplitude of the fundamental (pi) is determined in the same 
manner as in Fig. 2. The ratio of the amplitude of the harmonic (solid lines) to the square of the fundamental (dots) yields a 
value of 2F3go = 1.0(3) from Eqn.(2] from which we obtain r3go = 0.50(15). c, Dingle analysis, performed on a linear-log plot 
of the fundamental oscillations in contactless resistivity (denoted by solid circles, determined from the oscillation maxima and 
minima after dividing out Rt,i) at = 38° up to fields of 85 T, plotted versus 1/B. The 1/B = intercept obtained by fitting 
a straight line (purple) to this plot yields po = (1.6(5) x 10°) Hz. 



tiplicity of frequency components observed being a con- 
sequence of warping, bilayer splitting along with the ef- 
fects of magnetic breakdown that can produce extended 
orbits, and also potentially open orbits. This pocket 
is likely located at the nodal region, since photoemis- 
sion [l4 - 16 1 and high field heat capacity [TtJ experiments 
find a concentration in the density of states at the nodes 
of the superconducting wavefunction. Given that within 
translational symmetry breaking models proposed thus 



far (e.g. [111. Il2j]). the nodal pockets comprise hole carri- 



ers, we are faced with a dilemma as to how to explain the 
high field quantum oscillations with a negative Hall and 
Seebeck effect (3,113] m terms of these models, unless the 
Hall and Seebeck effects do not reflect the sign of car- 
rier in the Fermi surface pockets, or unless a significant 
change in electronic structure is induced by a magnetic 
field. Possible resolutions to this problem involve under- 
standing how a negative Hall and Seebeck effect can arise 
from a nodal pocket 27 1. 
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